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An approach to computing, withing the framework of distribution theory, the distributional val- 
ued energy-momentum tensor for the Schwarzschild spacetime is disscused. This approach avoids the 
"j^" problems associated with the regularization of singularities in the curvature tensors and shares COITI- 

ON ■ mon features with the by now standard treatment of discontinuities in General Relativity. Finally, 

' the Reissner-Nordstrom spacetime is also considered using the same approach. 

PACS numbers: 0420, 0470B, 0250N, 0420C,9760L 
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I. INTRODUCTION 

The energy-momentum tensor T associated to a point particle has an unambiguos character either if the particle 
is considered in a flat minkowskian spacetime or as a test particle in a curve background geometry. Obtained by a 
£SJ ' constructive procedure, invoking the minimal coupling principle or defined as the variational derivative of the point 
t"***" . particle action with respect to the metric, the result in an obvious notation, is 

o ■ 

T ^ = m 6 -^Mu^ (1) 
\F^9 dt 



If the particle is at rest with respect to the chosen reference frame, the only surviving component is the t—t component, 
O i.e., the one representing energy density, as it should be for a particle with zero momentum. 

The situation changes radically in general relativity if the point particle is acting as the gravitational field source. 
. The singular behaviour of the metric and related geometrical quantities at the particle location turns the mathematical 
meaning of the energy tensor rather imprecise. Nevertheless, on physical grounds or suggested by application of the 
variational principle leading to Einstein field equations, it seems plausible that the only non vanishing component 
■ of the energy- momentum tensor for a static point source is again the t—t component, with support at the particle 
K> \ position, and that is what it should be expected, whatever techniques are employed to find T. 

In the classical theory of gravitation one is led to consider the Einstein field equations which are, in general, 
quasilinear partial differential equations involving second order derivatives for the metric tensor. Hence, continuity of 
the first fundamental form is expected and at most, discontinuities in the second fundamental form, the coordinate 
independent statements appropriate to consider 3-surfaces of discontinuity in the spacetime manifolfd of General 
Relativity Q. 

The use of classical distribution theory, in order to incorporate singular information into General Relativity, was 
advocated some time ago However, it has been shown that the singular parts of the Riemann tensor must be 

of support on a submanifold of dimension d > 3, in order to have well defined tensor distributions for it and its 
contractions Q. In spite of this, several attemps have been made to regularize the Schwarzschild [j5],^| and Kerr- 
Newman spacetimes to obtain, after suitable limiting processes, the distributional energy-momentum tensors 
which play the role of point sources for these geometries with results that seem to depend on the regularization 
prescription ||. 

Besides the fact that the notion of a point source may has no precise mathematical meaning withing the framework 
of the nonlinear Einstein's theory of general relativity, some failures of previous works can be, to some extend, clearly 
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stablished and then circumvented. In order to reveal some of the weaknesses of the above forementioned regularization 
approaches, we will rederive the results of Ref. || following a different approach. Consider the line element 

ds 2 = h{r)dt 2 - h(r)- l dr 2 + r 2 d0 2 + r 2 sin 2 6 dcf 2 , (2) 

where 

9m 

h(r) = -l + — 6(r-e), (3) 
r 

with r < 2m, Q(r — e) is the Heaviside function and the limit e — > is understood. Equation with h as given in (|^) 
can be considered as a regularized version of the Schwarzschild line element in curvature coordinates. From equation 
(H), the calculation of the Einstein tensor proceeds in a straighforward manner. It gives 

G* t = G r r = --f-\(l + h), (4) 
r dr 

and 



, _ 1 d' 2 h 1 dh 



Next, taking the derivatives in the sense of distributions one finds 



G\ = G\ = ~2m S( - T C \ (6) 



and 



8(r — e) m d S(r - e) me d 

G e = G\ = -m 5 - —l-K* - £ ) = m 2 - —^ S ( r ~ £ )- ( ? ) 

e or r z dr 



We get in the limit e — > 



G ( t = G' r = -iG 9 e = -\G\ = -2m^M, (8) 

which is exactly the result obtained in Ref. || using smoothed versions of the Heaviside function S(r — e). This should 
be contrasted with what is the expected result 

G\ = -8nm5 3 {x)S%5° b . (9) 

Note that the strategy followed here to obtain equations (§|,0) suggest the possibility of attain a distributional meaning 
to curvature tensors by taking derivatives of the metric in the sense of distributions, an idea we shall pursue later. A 
second approach to the problem Q|, using the Kerr-Schild ansatz for the metric tensor, has been claimed to give the 
correct result (^), assuming that even under the regularization procedure the metric maintains its Kerr-Schild form. 
Nevertheless, it can be shown that in Kerr-Schild coordinates, without the constraint imposed by the Kerr-Schild 
ansatz, the same G a b as given by (^) is obtained. 

Several regularization schemes of this kind can be questioned on formal grounds. The unpleasant features are 
traceable, at least in part, to the fact that the metric obtained from the line element (|2|,[|) does not satisfies reasonable 
continuity requirements on the 3-surface r = e. Furthermore, as is well known for the Schwarzschild solution in 
curvature coordinates, r changes in character from a spacelike coordinate for r > 2m to a timelike one for r < 2m. 
Hence for values of r < 2m it is necessary to use a nonstatic system of coordinates. The discontinuity of the 
metric tensor and the change of the signature cast doubt over the previous results obtained from regularizing the 
Schwarzschild metric in these coordinates. 

In order to circumvent these problems we will consider a different approach which is closer to the classical theory 
of distributions view and that has a parallel into the theory of surface layers in General Relativity jlj . Assuming that 
the Schwarzschild geometry can be considered as being generated by an energy-momentum tensor with support in the 
singularity, one method of attacking the problem is to replace the point source by a uniform simple layer on a sphere 
of radius pQ. As po — ► the configuration tends to a point but for radius > po we always have an exact solution of 
Einstein's equations that looks exactly like the Schwarzschild solution. In a certain sense, we shall have to deal with 
a sequence of distributions, or with distributions depending on an arbitrary parameter po 
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The paper is organized as follows: in section II we briefly recall the construction of surface stress energy tensors for 
spherical thin shells, following the treatment of junction conditions of Israel Q. Some features of this construction are 
also discussed. In section III, using distributional techniques, the energy-momentum tensor on the whole manifold for 
a spherical thin shell is derived and showed to be in complete agreement with the results of previous section. Section 
IV is devoted to the point particle limit of the spherical shell, viewed as the limit of a sequence of distributions. 
Following the same approach, it is shown in section V, how these results are extended to the Reissner-Nordstrom 
spacetime. Finally, some concluding remarks are given in Section VI. 

We shall use geometrized units in which c~G~ 1, and the signature of the metric is (—,+,+,+) 



II. THE SURFACE LAYER FORMALISM. 



This section will be intented to apply the formalism of singular hypcrsurfaccs in order to construct the surface 
stress tensor of a spherical thin shell. We shall follow Lake notation ||, with minor differences in conventions. 

Let us consider a spherical thin shell E, a singular time- like hypersurface which divides the spacetime in two 
regions: the interior region V~, described by flat minkowskian geometry and V + , the exterior spacetime, chosen as 
a spherically symmetric solution to Einstein held equations. Both regions V r± will be decribed by a single set of 
isotropic coordinates, x a = (£, p, 9, 4>) with respect to which the line element is 

ds 2 = -A 2 dt 2 + B 2 {dp 2 + p 2 d6 2 + p 2 sin 2 6 dcf 2 ) (10) 

where for the exterior geometry A = A(p) and B = B(p), are known functions of the radial isotropic coordinate and 
for the interior region, A = A(po) = Aq and B — B(po) = Bq are constants. The surface layer is represented in the 
chosen frame by 

f(p) = p-po = (11) 

therefore, the components of the metric tensor matches naturally and smoothly across the shell, corresponding to the 
well defined character of the intrinsic geometry of E, in fact, the induced metric on E takes the form 

ds| = ~A 2 dt 2 + B' 2 pl(d0 2 + sin 2 9 d(j) 2 ) (12) 

where intrinsic coordinates = (t, 9, <f>) are used and the 3-metric elements induced are 

0x a dx^ 

For future purpose we note here that the relation between coordinates of and intrinsic coordinates on E are 

dr a 

%-t ("> 

On the other hand, the extrinsec geometry of E is not well defined, i.e., the extrinsec curvature tensor, K = — ^£ n g 
induced on the layer is diferent as evaluated from its embedding in V + or V~ and the surface stress energy tensor is 
related to the discontinuity of K through the Lanczos equation 

8nS^ = [K^] - SjTr [K] i, j = (i, 9, <f) (15) 







where square brakets denotes a discontinuity across de layer as evaluated in V + and V~ , i.e., [/] = / + — f~. Our 
task is to perform the calculations to obtain the surface stress energy tensor of the layer. From its definition it is easy 
to see that 



dx a dx 



(16) 



n a being the components of n, the outward unit normal field to the 3-surface E. Using (11) and (10) the normal is 
found to be 
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BS" a (17) 



With the usual formula for the covariant derivative and equations (|1J) and (|17|), equation ( |16[) yields the simple 
relation 

K tj = BY% , (18) 

where Y stand for the usual Christoffel symbols. Thus, raising indexes, and performing the calculations for V + and 
for V~ we obtain 

[K\] = ~, [K%] = [K%] = - * (19) 

where and hereafter primes denotes differentiation with respect to p and it is understood that quantities must be 
evaluated at E after performing the derivations. Finally, using equation ([if]) for the surface energy tensor, we arrive 
to the simple results 

s '< = S5IT (20) 

and 

ce o<j> _ 1 ( AB )' (21 ) 

Let us apply the above results to obtain the surface stress tensor for a spherical layer of total gravitational mass 
m. Therefore A and B will be the metric elements corresponding to Schwarzschild vacumm solution in isotropic 
coordinates 

As = (l-^)(l + ^ p )- 1 in (22) 



and 



B S = (1 + ^) 2 in V+ (23) 
2p 



where as usual, m is the gravitational mass as measured at infinite. From (20) and (|2l|), S can be readily obtained as 



S\= -f- 2 {l + ^T* (24) 



and 



v ISnpl 2p 2p Q 

The surface energy density of the layer tr, defined as a = SijU l u J , where u % are the components of the intrinsic tangent 
vector, u 1 — A~ 1 Sl, is given by a — S\. The energy of the shell can be evaluated by integrating a over the proper 
area of the 2-sphere t= const, and p = const. A simple calculation gives 



Sau l u 3 da = / a da — m H . (26) 

1 2p 



III. THE DISTRIBUTIONAL APPROACH TO THIN SHELLS. 

Consider again the line clement 

ds 2 = -A 2 (p)dt 2 + B 2 (p)(dp 2 + p 2 d9 2 + p 2 sin6 dc/) 2 ), (27) 
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A{ P ) 



{ 2po> 



) 



(p < Po) 
{P > Po) 



(28) 



B(p) 




(P < Po) 
(P > Po) 



(29) 



which corresponds to the geometry produced by a very thin spherical shell of isotropic radius po'- inside the sphere, the 
metric is equivalent to the flat-space Minkowski metric and outside it agrees with the Schwarzschild solution written 
in isotropic coordinates, provided we restrict ourselves to values of po satisfying the condition p > \vn. 



From 



we obtain for the Einstein tensor 



G p , 



2B"Bp- B' 2 p + 4B'B 
B' 2 A p + 2B'AB + 2A'B'B p + 2A'B 2 



AB 4 p 

B'AB + B" ABp + A'B' 2 - B' 2 Ap + A"B 2 p 



G% = G\ = 



AB 4 p 



(30) 
(31) 
(32) 



Because of the linear character of the theory of distributions, the differential equations satisfied by them should be 
linear. The above expressions can not be written in the form of linear differential operators acting on A and B nor 
even on some funtions of these. To what extent can the Einstein tensor be regarded as a distribution? What we shall 
consider here is an alternative that one should expect to find useful in many nonlinear problems. In classical analysis, 
A and B are locally integrable infinitely diffcrentiable functions of variable p except at p = po where A' and B' are 
discontinues and A" and B" fail to exist. However, differentiating A and B in the sense of distributions, A' and B' 
are now distributions with jump discontinuities at £ with resulting fS(p — pp) terms in the expressions for A" and 
B" , with / the corresponding jump of the first derivative. From equations (|30-32), this procedure provides us with a 
physically sensible well defined distributional tensor G a b given by 



1 



p 2 AB- 



-2m(\ 



2 Pl 



-)S(p-po), 



G" p = 0, 
G e B = Gt = 



2p 2 AB 3 po 



6{p-po), 



(33) 
(34) 

(35) 



with support on the 3-surface S. This confirm the interpretation of equations (|27[|29j) as the geometry whose source 
corresponds to a uniform surface layer density spread on the 3-surface S, i.e., to an energy-momentum tensor T with 
support on E: 



K = o, 



1 m . 

.___ 2m (l-—)5(p-p ), 



rp<t 



16irp 2 AB 3 po 



Hp - po), 



(36) 
(37) 

(38) 



where use have been made of the Einstein equation G = 87rT. We recall that, the linear character of the theory 
of distributions avoids a more rigorous justification of the above result. However, its correctness follows directly 
from the relation between the surface energy- momentum tensor S a b on the 3-surface S, equations (^J,H|), and the 
energy-momentum tensor T a b on the whole manifold ([3(||3^) , namely Jl^] 



S a b = lim 

£^0 



T\dn 



(39) 



with dn — Bdp being the proper distance element measured perpendicularly through the 3-surface E. 
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The energy-momentum tensor T can be cast in the form of the energy-momentum tensor of an anisotropic fluid 

T% = Qu a u h + P 1 _h\ (40) 

where £ is the proper energy density, u a are the components of the unit timelike four velocity, P±_ is the tangential 
pressure and h a b is the projection tensor onto the subspace orthogonal to u and n, 

h\ = 81 + u a u b - n a n b (41) 

with 

(I-^)HP-Po) (42) 



s 4irp 2 AB 3 y 2p 
and 

P - = ^AB^ p -^ (43) 
Note that the calculation of the energy of the layer as measured by an observer with velocity u gives 

£ = J T ab u a u b d 3 a = J (d 3 a = m(l + ^-), (44) 

where d 3 a is the volumen element of the rest space of the observer, in complete agreement with (p(f). Furthermore, 
the Tolman formula for the total energy of a static and asymptotically flat spacetime 

£t = J (T p p + T% + T\ - T\)^gd 3 x, (45) 

with g the determinant of the four dimensional metric and d 3 x the coordinate volume element, gives 

E T = m, (46) 

as it should be. Interesting enough, note that the difference between the total energy as measured in infinite and the 
energy of the shell, is — m 2 /2po, which is the expression for the gravitational newtonian energy of a thin shell; recall 
however that the radial coordinate po is not the proper radius of the shell. 

IV. THE POINT PARTICLE LIMIT AND THE SCHWARZSCHILD GEOMETRY. 

Now we consider the limit in which the spherical layer tends to a point particle and let us comment on the limiting 
procedure. It should be noted that the value of m/2po must be keep fixed and satisfying the condition m/2po < 1 
in order to have a well defined metric tensor, equations ( p7| , p8 29) for all p. Furthermore, the static character of 



this spacetime reduce the distributional evaluation of lim Po _->o G a b to a three-dimensional problem on the 3-surfaces 
t = const. For m/2po fixed, it turns out that the sequence of distributions G a b is such that lim Po ^o (G a b , $) exist for 
each $. Let us ilustrate the evaluation of linip^o G a b with the explicit calculation of lim po _>o G t where the condition 
m/2po — fixed is understood. We have 

(G\, *(*)>, = (-^2to(1 - g)8(p - Po)M*)) pM 

r°° fx r 2 * i m 

= -/ dp d9 d<f>p 2 B 3 sm6 __2m(l- —)S(p-p Mp,9 ,</>) 
Jo Jo Jo p 1 AB i 2p Q 

= -2m(l+— ) / del d(j)sm9 ${po 7 9,<f>). 
2po Jo Jo 



Taking into account that 



,. . to / , m 

lim m(l + — ) = lim to + 2p Q ( — )■ 
Po^o 2p Po^o \ 2p 



G 



it follows that 

lim (G* f , $(£))., = lim -2m(l 



Po^O 



lim 



7T /•2-7T 



d<9 / d^sind [$(0) + ($(p,6,<f>) - $(0)] 



-87rm(l + )$(0)-2m(l + — ) / d0 / . 

2po 2p Jo Jo 



>sin6> (<f>(p,6,<p) ~ $(0) 



-87rm$(0) = (-8tttoc5 3 (£),$(£))_ 



Hence 



lim G\ = -8irm5 3 {x) 



(47) 



The evaluation of the limit for the remaining non zero components of G a b proceeds in the same manner. The results 
are 



lim G% = lim G\ = 

It them follows that this limit is a well defined distribution given by 

G\ = -8irm5 3 (x)6%6%, 

which implies via Einstein equations, 

T\ = -m6 3 (x)6S5° b , 

as expected. 



(48) 

(49) 
(50) 



V. THE REISSNER-NORDSTROM SPACETIME. 

Next, let us consider the case in which the shell has a uniformly distributed charge e. Hence, the V + spacetime line 
element is given by the Reissner-Nordstrom exterior solution to the Einstein-Maxwell equations, appropiately written 
in isotropic coordinates, i.e., 



Arn — 



1 2p> 4p 2 



V+ 



(51) 



and 



It follows from (§9) and (§|) that 



and 



~ / fit/ x o c 



Sirpt 



2m(l + — )-- 
2p Po 



b rn(Po) 



2 2 

to - e 



b rn(po) a rn(po)- 



(52) 



(53) 



(54) 



Some comments are in order. First observe that in the extreme Reissner-Nordstrom solution, e = ± to, the tangential 
stresses vanish. This is explained if we realize that e = ± to is precisely the Bonnor condition []l3f for a selfgravitating 
pressureless charged fluid to be in equilibrium, therefore in this case we are dealing with a dust thin shell. Secondly, 
note that the calculation of the energy contained in the charged layer, integrating (^3|) as was done in (^6|), gives 
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£ = m(l 



m , 
Wo 



e 



(55) 



the last term in this equation representing the (negative) electrostatic contribution to the energy. Finally, in the limit 
e = these equations reduce to the neutral expected. 

Following the same approach of section III, the energy-momentum tensor for the spherical thin shell with a uniformly 
distributed charge e can be obtained. It is given by 



T , 



1 



8irp 2 A RN B 3 



RN 



0_ 

(1- 



m \ _|_ e 



2p < 



4PS 



2 1 



2m (1 



e 
Po 



S(p- Po) 



8lip 4 B 4 



-6(p- po), 



RN 



rpp 
P 



8irp 4 B 4 



-&(p- po), 



UN 



n = t\ = 



+ 



16np 2 A RN B RN 



A) 



-$(p-Po) 



8irp i B 4 



-Q(p-Po). 



(56) 
(57) 

(58) 



RN 



From these results one can verify, using equation (|39j) , that the energy-momentum tensor S on the hipersurface S, 
equations (53 54|) is in fact the "delta- like singularity" at £ [[l2| of the energy- momentum tensor T given by (|5^-p8|). 
Note that the non-delta- like electromagnetic part of (56-58|) is traceless. Furthermore, from ( fi5"l) and ( p6[j5^ ), it follows 
the remarkable result £t = m, i.e., £t no contains electromagnetic contributions. 

The limiting case of a point particle runs over the same lines given in Section IV. The result is 



-mS 3 (x) 



8np 4 B 4 



■e(p), 



RN 



rpp 
P 



8np 4 B 4 



-e(p), 



T U = rp<f 



RN 

e 2 

8np 4 B 4 



-e(p). 



(59) 



RN 



as expected. 



VI. CONCLUSIONS AND REMARKS. 



We have shown that a succesfull approach for dealing with curvature tensor valued distribution is to first impose 
admisible continuity conditions on the metric tensor, and then take its derivatives in the sense of classical distributions. 
The distributional meaning is then equivalent to the junction condition formalism. Afterwards, through appropiatc 
limiting procedures, it is then possible to obtain well behaved distributional tensors with support on submanifolds 
of d < 3, as we have shown for the energy- momentum tensors associated with the Schwarzschild and Reissner- 
Nordstrom spacetimes. The above procedure provides us with what is expected on physical grounds. However, it 
should be mentioned that the use of Colombeau's new generalized functions p[ in order to obtain distributional 
curvatures |jl^,^6|, may renders a more rigorous setting for discussing situations like the ones considered in this paper. 
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